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Abstract 

An Ising model with ferromagnetic nearest-neighbor interactions Ji (Ji > 0) and random 
next-nearest- neighbor interactions [+J2 with probabihtyp and — J2 with probabiHty (1— p); 
J2 > 0] is studied within the framework of an effective-field theory based on the differential- 
operator technique. The order parameters are calculated, considering finite clusters with 
n = I, 2, and 4 spins, using the standard approximation of neglecting correlations. A phase 
diagram is obtained in the plane temperature versus p, for the particular case Ji — J2, 
showing both superantiferromagnetic (low p) and ferromagnetic (higher values of p) orderings 
at low temperatures. 
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Many magnetic compounds are well-described in terms of theoretical models characterized 
by a competition of nearest-neighbor and next-nearest-neighbor interactions. As typical exam- 
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pies, one has EuxSri_xS [l|, y] and FexZni„xF2 [3|, which may present various low-temperature 
magnetic orderings, depending on its parameters, like the strength of these interactions and the 
concentration of magnetic ions x. The simplest model for such systems consists of an Ising model 
with competing uniform interactions, ferromagnetic Ji (or antiferromagnetic — Ji) (Ji > 0) on 
nearest-neighbor, and antiferromagnetic — J2 (J2 > 0) on next-nearest-neighbor spins (to be 
denoted hereafter as Ji — J2 Ising model). 

The Ji — J2 Ising model on a square lattice has attracted the attention of many workers, b eing 



investi gate d through several approximation methods 
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23]. The various possible phases are the paramagnetic 



(P), for high temperatures, whereas for low temperatures one may have the ferromagnetic (F), 
antiferromagnetic (AF) and superantiferromagnetic (SAF) phases; the later is characterized 
by alternate ferromagnetic rows (or columns) of oppositely oriented spins. In the absence of a 
magnetic field, it can be shown that one has a symmetry with respect to the sign of the nearest- 
neighbor interactions, i.e., the ferromagnetic and antiferromagnetic states are equivalent Q]. 
At zero temperature (ground state), one has two ordered states depending on the value of the 
frustration parameter a = {J2/J1), namely, the F (0 < a < 1/2) (in an equivalent way, an AF 
state occurs instead of the F one, if the sign of nearest-neighbor interactions is changed) and the 
SAF (a > 1/2) phases. On the other hand, its phase diagram for finite temperatures has been 
the object of some controversies. A continuous critical frontier between the P and F phases, for 
< a < 1/2, with the critical temperature Tc{a) decreasing and approaching zero for a — > 1/2, 
is well-accepted nowadays. However, some characteristics of this phase diagram for a > 1/2 are 
polemic; in particular, the presence of a first-order transition between the P and SAF phases, 
as well as a tricritical point characterized by the coordinates {at, Tt), beyond which a continuous 
critical frontier occurs, has been the object of some debate. Most of the works indicate that 
within the range 1/2 < a < at one should get a first-order phase transition between the P and 
SAF phases, whereas for a > at this transition should become continuous. Some studies using 

12, Q, lil claim a first-order line in the range 0.5 < a < 1.2, i.e.. 



cluster-variation methods 



at = 1.2. However, recent results of Monte Carlo simulations [241 . 



partition function 



26l |. analysis of zeros of the 



25l |. and effective-field theory (EFT) [27|] suggest that at < 1, in such a way 



that a continuous transition between the P and SAF phases should occur for a = 1. 
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In the present work we introduce a disordered version of the Ji — J2 Ising model on a 
square lattice. It is important to mention that due to the symmetry with respect to the sign 
of the nearest-neighbor interactions, for which a change of sign takes F ^ AF, a disorder that 
produces a simple change of sign in Ji is not expected to bring any novel physical properties 
5]. A disordered version of the Ji — J2 Ising model has already been defined in Refs. [3, 13 
by considering the same bond disorder, P{Jij) = {l/2)[6{Jij — Ji) + 5{Jij — J2)], for both 
nearest- and next-nearest-neighbor pairs of spins, which were introduced with the same sign 
in the corresponding Hamiltonian. In this analysis, the authors have investigated particular 
choices for the ratio J2/J1, through numerical simulations of the model on a square lattice, 
and have found as its main results: (i) An interesting saturation, with respect to the size of 
the system, of the specific-heat maxima; (ii) Critical exponents following the weak-universality 
scenario, according to which critical exponents change with the disorder in a certain way to 
keep some ratios (e.g., (3/u and ^/v) unchanged at their corresponding pure values, similarly to 
what happens for other disordered magnetic models. Herein, we introduce a different disordered 
version of the Ji — J2 Ising model, by defining a disorder only in the next-nearest-neighbor 
interactions, in such a way that we shall consider the Hamiltonian, 



where the summations are over nearest-neighbor (nn) and next-nearest-neighbor (nnn) pairs of 
spins, with ferromagnetic ( Ji > 0) and random J- interactions, respectively. The later follow 
a bimodal probability distribution. 



V{jlf)=p5{j\f -h) + {l-p)5{jlf +,h) (J2>0). (2) 



The particular case p = of this model corresponds to the pure Ji — J2 Ising mode 
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28( 1 . The introduction of the above type of 



disorder allows for the possibility of a richer critical behavior, characterized by a competition 
between different types of orderings as one varies p, for a fixed ratio of the interaction strengths. 
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a = {J2/ Ji)- In addition to that, from the point of view of physical reahzations, the parameter 
p, which is usuahy associated with the concentration of a given chemical element, is much easier 
to be varied experimentally than the ratio of interaction strengths, a. In the present analysis 
we restrict ourselves to the case J2 = Ji, i-e., a = 1. 

Herein, we treat this model by using the EFT method, which has been shown recently to be 
very useful in the pure case (i.e., p = 0) [27]. The EFT technique is based on rigorous correla- 
tion identities, as a starting point, and applies the differential-operator technique developed by 
Honmura and Kaneyoshi [31]. Within this method, the averages of a general function A{{n}), 
involving spin variables, is obtained by 



usual canonical thermal average over the surrounding of the cluster. In the case n = 2 (see 
Fig. [l]) the trace Tr{„| is calculated over spins Si and ^2, whereas (• ■ •) is to be considered over 
C7i,(T2, • • • ,aio. 

We have applied the EFT method for clusters of with n = 1,2, and 4 spins (to be denoted 
hereafter EFT-n) and below we describe the case of EFT-2, illustrated in Fig. [TJ The Hamilto- 
nian for this cluster is given by 




(3) 



where the partial trace Tr^^j is calculated exactly over the set {n} of spin variables that belong 
to the finite cluster specified by the multisite spin Hamiltonian 7^{„} and (• • •) indicates the 



n2 



JiSi ■ 5*2 + Siai + 5202 , 



(4) 



with 
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(5) 



i=l 



j=4,5,7,8 



and 
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i=2,3,9,10 
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We have divided the square lattice into two sublattices A and B, each of them defined by 
alternating lines (or columns) in such a way that in Fig.[T]one has niA = {{Si))j(2) = {{S2))j(2) = 
((cTi))j{2) (i = 1,6) andm^ = ((c"i))j{2) {i = 2,3,4,5,7,8,9,10), where (••)j(2) denotes an average 
over the disorder; in the F (SAF) phase one has that ttia = ttib = rn {niA = —ruB = m). At 
zero temperature [m = 1), the ground state of the Hamiltonian of Eq. ([T]) is ferromagnetic for 
p = 1 [see Fig. [D^a)], whereas for p = and a = 1 one gets the superantiferromagnetic state 
[shown in Fig. [IJb)]. In each case, a continuous phase transition occurs at a finite temperature 
between the corresponding low-temperature state and the paramagnetic state. To the best of 
our knowledge, theoretical works to investigate the critical behavior of this model for < p < 1 
have never been carried and this represents the purpose of the present work. 

The parameter m may be obtained by evaluating the average magnetization per spin in 
sublattice A, itia = {{^{Si + 5*2))) j(2), calculating the inner trace in Eq. ^ over spins Si, S2 = 
±1, which for the Hamiltonian of Eq. dH) yields, 



m = {{ ^ sinh(ai + a2) ^ ^ \\ ^ .7) 

cosh(ai + 02) + exp(-2Ki) cosh(ai - 02) / / j(2) 



where oi = — /3ai, 02 = — /3a2, and Ki = (5Ji. Using the identity exp(aZ)a; + bDy)g{x,y) 

d_ 



g{x + a,y + b), where — '§71 (/^ = 2;, y) is the differential operator, Eq. ([7]) becomes 



m = {{eyiv{aiD^ + a2Dy)))j^2) g{x,y)\^y^Q , (8) 



with 



I X sinh(x + y) 

cosh(x + y) + exp(— zivi) cosh(x — y) 



Applying the van der Waerden identity for the exponentials containing Ising variables, i.e., 
exp(acrj) = cosh(a) +(Tj sinh(a) (cjj = ±1), the right-hand side of Eq. ([8]) can be written exactly 
in terms of multiple spin correlation functions; however, it is clear that if one tries to treat 
exactly all correlation functions, the problem becomes unmanageable. In this work, we use a 
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decoupling procedure that ignores all high-order correlations on the right-hand side of Eq. ([S]), 
i.e., 

((cJi -aj-'-ar C7p))j(2) ^ {{(^i)) j{2) ■ {{(Tj))j(2) ■ ■ ■ ((o-/))j(2) • {{crp))j(2) , (10) 

where i^jy^---^ly^P- The present procedure neglects correlations between different spins, 
but takes into account relations such as (di)^ = 1, whereas in the usual mean- field approximation 
both self- and multi-spin correlation functions are neglected. Applying the approximation of 
Eq. (fTOll in Eq. ([8]), one gets the equations of state. 



Y,Ag\,iT,p)m''~^\ (11) 



m = 

k=0 



where, due to the symmetry of the Hamiltonian, the even powers of m do not occur. In the 

(V) 

equation above, the coefficients A2i^^^{T,p) depend on the boundary conditions shown in Fig.[Tl 
through the index F = F, SAF; in the present EFT-2 analysis one has that the highest power 
occurring on the right-hand side is defined by /c = 4. In general, the coefficients 712^+1(^5?*) 
depend on the size n of the cluster used in EFT-n, getting more complicated for larger clusters; 
in addition to that, k increases with n, e.g., k = 3 for n = 1 and k = 5 for n = 4, in such a way 
that the procedure may get very tedious for large clusters. 

The borders of continuous phase transitions between the ordered (F and SAF) phases and 
the high-temperature disordered (P) phase may be obtained by considering the limit m ^ in 
Eq. (fTT|) . which leads to 

APiT,,p) = l. (12) 

By solving Eq. (|12|) numerically, one finds Tc{p) for both phases F and SAF, as shown in the 
phase diagram of Fig. [2j For p = one has a continuous transition between the SAF and P 
phases at {ksTc^O)/ Ji) = 2.263, which is in g^ood agreement with the value (A;bTc(0)/Ji) = 2.083 
estimated from Monte Carlo simulations 



10|,l2|, as well as from a computation of the zeros of 



7 



the partition function 25|]. The critical temperature associated with the SAF phase decreases as 
p increases, in such a way that one gets Tc{pci) = 0, for pd = 0.075. Due to strong frustration 
effects, the system does not present long-range order in the interval pd < p < Pc2, where 
Pc2 — 0.471. For p > pc2 one gets that the critical temperature associated with the F phase 
increases monotonically, leading to (A;bTc(1)/Ji) = 6.937. Besides the results calculated within 
EFT-2 shown in Fig. [21 we have also computed these critical frontiers within EFT-1 and EFT-4. 
Only in the EFT-1 case, the critical frontier SAF-P comes as a first-order transition; we believe 
this to be a spurious result (attributed to the smallness of the cluster), since it is in contrast 
with those obtained from EFT-2 and EFT-4 (which yield continuous phase transitions), as well 
as with those of Monte Carlo simulations for the particular case p = From EFT-4 

one gets that (A;bTc(0)/Ji) = 2.310 and pd = 0.088, indicating an enlargement of the SAF 
phase, whereas for the F phase one has Pc2 — 0.474 and {kBTc{l)/ Ji) = 6.750. It should be 
mentioned that our estimates for (A;b?c(1)/Ji), from EFT-1, EFT-2, and EFT-4, suggest a 
slow convergence towards those from low- and high-temperature series expansions, which yield 
(A;bTc(1)/ Ji) = 5.260 and (ksTcil) / Ji) = 5.257, respectively [34], or to the more recent one. 



221. 



obtained from an analytical expression for the interface free energy, (kBTc{l)/ Ji) = 5.376 
In a similar way, our results suggest a decrease in the gap {pc2 — Pel), although it is not possible 
to conclude whether the phases SAF and F should meet at zero temperature. 

It is important to mention that the EFT method, which replaces averages over products of 
spins by the respective products of their averages, neglects correlations and, as a consequence, 
the associated critical exponents are mean-field-like. Therefore, this method is not suitable 
for investigations of possible changes in the critical exponents with respect to variations of 
important parameters of the problem (like R and p), which represents a major question in 
its p = particular case, the Ji — J2 Ising model. Furthermore, since the EFT-n procedure 
consists in treating a cluster of n spins exactly, whereas the interaction of this cluster with 
its surrounding are treated as an average (mean- field-like) , the method becomes unmanageable 
for increasing values of n. However, this technique presents the advantage of providing phase 
diagrams that may be qualitatively correct as a whole, and in some cases it may yield critical 
points with a good accuracy, in spite of considering clusters with small n, like in the following 
examples: (i) The tricritical point in the phase diagram of the Ji — J2 Ising model on the 



square lattice, which was estimated to occur for {J2/JxL< 1 within EFT [27l], is in agreement 



with other methods, like Monte Carlo simulations 



24 



26l | and analysis of zeros of the partition 



function [25]; (ii) The present estimate of the critical temperature for p = 0, (A;bT'c(0)/Ji) 



2.263, yields a relative discrepancy of typically J 
2.083, estimated from Monte Carlo simulations 



^ with respect to the value (/cbT'c(O)/ J; 
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24| and zeros of the partition function 



25(. In the present problem, a full phase diagram was obtained for the case (J2/J1) = 1, 
where an elimination of the SAF phase occurs as one increases the concentration of second- 
neighbor ferromagnetic bonds. Therefore, according to the present EFT approach, only for 
low values of p one gets a sufficient concentration of second-neighbor antiferromagnetic bonds 
in order to yield the SAF ordering at low temperatures. As the concentration p increases, 
the competition between ferromagnetic couplings (both from the nearest-neighbor pairs and 
second-neighbor pairs with probability p) and the second-neighbor antiferromagnetic ones [with 
probability (1 — p)] destroys the SAF ordering at the concentration pd- Above this value, the 
competition between ferro- and antiferromagnetic interactions leads to a sufficient amount of 
frustrations that could, in principle, favor a type of spin-glass ordering at low temperatures. 
However, similarly to many other frustrated two-dimensional systems, which do not exhibit a 
spin-glass phase for finite temperatures, herein the paramagnetic phase dominates up to zero 
temperature, in the interval pd < p < Pc2- For p > pc2, one has sufficient ferromagnetic bonds 
leading to the F ordering. Obviously, the present estimates of pci and pc2 are expected to change 
under other approximation techniques (e.g., Monte Carlo simulations) and one cannot rule out 
the possibility of a critical frontier separating the phases SAF and F at finite temperatures. 

To conclude, we have introduced an Ising model characterized by ferromagnetic nearest- 
neighbor interactions Ji and random ibJ2 next-nearest-neighbor interactions [+J2 with prob- 
ability p and —J2 with probability (1 — p)]- This model, which generalizes the Ji — J2 Ising 
model (herein corresponding to the case p = 0) is expected to be relevant for some diluted mag- 
netic compounds, characterized by a competition between nearest- and next-nearest-neighbor 
interactions. We have studied the phase diagram of this model in the particular case Ji = J2 
within the framework of an effective-field theory based on the differential-operator technique, 
considering finite clusters with n = 1,2, and 4 spins. We have found a low-temperature sublat- 
tice ordering associated with a superantiferromagnetic phase for small p {0 < p < pd, where 
Pd = 0.088 in the case n = 4), as well as a ferromagnetic phase for larger values of p (p > Pc2, 
where Pc2 — 0.474 in the case n = 4), whereas in the interval pd < p < Pc2 the system does 
not present long-range order. Other approximation methods, such as renormalization-group and 
Monte Carlo simulations, should be used to obtain further information on the critical behavior 
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of this model. 
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Figure Captions 

Fig. 1: A cluster with n = 2 (represented by spins Si and S2) and its surrounding (spins 
cri,cJ2,--- ,cio)) for the model defined in Eq. ([1]), in its ferromagnetic (a) and superantiferro- 
magnetic (b) ground states. The full lines represent the interactions of spins Si and ^2, which are 
taken into account exactly in Eq. dSj) , whereas the dashed lines represent the usual square-lattice 
structure. 

Fig. 2: Phase diagram of the model defined in Eqs. ([T]) and ([2]) for J2 = Ji within EFT-2. 
The phases are the paramagnetic (P), ferromagnetic (F), and superantiferromagnetic (SAF), 
as described in the text. 
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FIG. 1: A cluster with n = 2 (represented by spins and 5*2) and its surrounding (spins cti, (72, • • • , fio), 
for the model defined in Eq. Ilj, in its ferromagnetic (a) and superantiferromagnetic (b) ground states. 
The full lines represent the interactions of spins 5*1 and 5*2, which are taken into account exactly in 
Eq. ([3]), whereas the dashed lines represent the usual square-lattice structure. 




FIG. 2: Phase diagram of the model defined in Eqs. ([T]) and ^ for J2 — Ji within EFT-2. The phases 
are the paramagnetic (P), ferromagnetic (F), and superantiferromagnetic (SAF), as described in the 
text. 



